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Abstract. We consider an ergodic Schrodinger operator with magnetic field 
within the non-interacting particle approximation. Justifying the linear re- 
sponse theory, a rigorous derivation of a Kubo formula for the electric con- 
ductivity tensor within this context can be found in a recent work of Bouclet, 
Germinet, Klein and Schenker [BoGKSj . If the Fermi level falls into a region of 
localization, the well-known Kubo-Streda formula for the quantum Hall con- 
ductivity at zero temperature is recovered. In this review we go along the lines 
of |BoGKSl but make a more systematic use of noncommutativc L^^-spaces, 
leading to a somewhat more transparent proof. 
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1. Introduction 

In [BoGKSj the authors consider an ergodic Schrodinger operator with magnetic 
field, and give a controlled derivation of a Kubo formula for the electric conductiv- 
ity tensor, validating the linear response theory within the noninteracting particle 
approximation. For an adiabatically switched electric field, they then recover the 
expected expression for the quantum Hall conductivity whenever the Fermi energy 
lies either in a region of localization of the reference Hamiltonian or in a gap of the 
spectrum. 
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The aim of this paper is to provide a pedestrian derivation of [BoGKSj 's result 
and to simphfy their "mathematical apparatus" by resorting more systematically to 
noncommutative L^-spaces. We also state results for more general time dependent 
electric fields, so that AC-conductivity is covered as well. That von Neumann 
algebra and noncommutative integration play an important role in the context of 
random Schrodinger operators involved with the quantum Hall effect goes back to 
Bellissard, e.g. |a iBESl [SBll [SB2] . 

The electric conductivity tensor is usually expressed in terms of a "Kubo for- 
mula," derived via formal linear response theory. In the context of disordered 
media, where Anderson localization is expected (or proved), the electric conduc- 
tivity has driven a lot of interest coming for several perspective. For time reversal 
systems and at zero temperature, the vanishing of the direct conductivity is a 
physically meaningful evidence of a localized regime |FS| lAGj . Another direction 
of interest is the connection between direct conductivity and the quantum Hall 
effect |ThKNN[ [Stl E KB iBESl EG] . On an other hand the behaviour 

of the alternative conductivity at small frequencies within the region of localiza- 
tion is dictated by the celebrated Mott formula jMD| (see |KLP| IKLM| IKMj for 
recent important developments). Connected to conductivities, current-current cor- 
relations functions have recently droven a lot of attention as well (see [BH, ,CGHj 
and references therein). 

During the past two decades a few papers managed to shed some light on these 
derivations from the mathematical point of view, e.g., [0 IKu) lB| INBi lAvSS) IBES[ 
[SBll [SB2l lAGl iNal lESl ICo.TMl ICoNP) . While a great amount of attention has been 
brought to the derivation from a Kubo formula of conductivities (in particular of 
the quantum Hall conductivity), and to the study of these conductivities, not much 
has been done concerning a controlled derivation of the linear response and the 
Kubo formula itself; only the recent papers [SB2l lESl IBoGKSI IGo.TMl IGoNPj deal 
with this question. 

In this note, the accent is put on the derivation of the linear response for which 
we shall present the main elements of proof, along the lines of jBoGKS] but using 
noncommutative integration. The required material is briefly presented or recalled 
from IBoGKSj . Further details and extended proofs will be found in [Do . We 
start by describing the noncommutative L^-spaces that are relevant in our context, 
and we state the properties that we shall need (Section [2]). In Section [3] we de- 
fine magnetic random Schrodinger operators and perturbations of these by time 
dependent electric fields, but in a suitable gauge where the electric field is given 
by a time-dependent vector potential. We review from iBoGKS, the main tools 
that enter the derivation of the linear response, in particular the time evolution 
propagators. In Section |4] we compute rigorously the linear response of the system 
forced by a time dependent electric field. We do it along the lines of jBoGKS] 
but within the framework of the noncommutative L^-spaces presented in Section [5] 
The derivation is achieved in several steps. First we set up the Liouville equation 
which describes the time evolution of the density matrix under the action of a time- 
dependent electric field (Theorem 14. ip . In a standard way, this evolution equation 
can be written as an integral equation, the so-called Duhamel formula. Second, we 
compute the net current per unit volume induced by the electric field and prove 
that it is differentiable with respect to the electric field at zero field. This yields in 
fair generality the desired Kubo formula for the electric conductivity tensor, for any 
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non zero adiabatic parameter (Theorem 14.61 and Corollarv l4.7p . The adiabatic hmit 
is then performed in order to compute the direct / ac conductivity at zero temper- 
ature (Theorem 14.81 CoroUary 14.91 and Remark 14. lip . In particular we recover the 
expected expression for the quantum HaU conductivity, the Kubo-Sti^eda formula, 
as in [Bl IBES| . At positive temperature, we note that, while the existence of the 
electric conductivity measure can easily be derived from that Kubo formula |KM) . 
proving that the conductivity itself, i.e. its density, exists and is finite remains out 
of reach. 

Acknowledgement. We thank warmly Vladimir Georgescu for enlightening dis- 
cussions on noncommutative integration, as well as A. Klein for useful comments. 

2. COVARIANT MEASURABLE OPERATORS AND NONCOMMUTATIVE LP-SPACES 

In this section we construct the noncommutative L^-spaces that are relevant for 
our analysis. We first recall the underlying Von Neumann alegbra of observables 
and we equip it with the so called "trace per unit volume" . We refer to [Dl ITe| for 
the material. We shall skip some details and proofs for which we also refer to [Do] . 

2.1. Observables. Let H be a separable Hilbert space (in our contexts = L2(M'^)). 
Let Z be an abelian locally compact group and U — {Ua}aez a unitary projective 
representation of Z on T-L, i.e. 

• UalJb = i{a,b)Ua+b, where ^(a, &) £ C, \£,{a,h)\ = 1; 

• Ue = l; 

Now we take a set of orthogonal projections on H , x := {Xa}a^z, Z — > B{'H). 
Such that if a / 6 XaXb = and J2aez Xa Furthermore one requires a 
covariance relation or a stability relation of x under U i.e UaXbU* = Xa+b- 

Next to this Hilbertian structure (representing the "physical" space) , we consider 
a probability space (fi, J", P) (representing the presence of the disorder) that is 
ergodic under the action of a group r = {Ta}aez, that is, 

• \/a G Z, Ta : — !• is a measure preserving isomorphism; 

• Va, 6 G Z, TaOTb^ Ta+t; 

• Te = 1 where e is the neutral element of Z and thus t^^ — r^a, E Z; 

• If ^ G is invariant under r, then P{A) = or 1. 

With these two structures at hand we define the Hilbert space 

n= nd¥{Lo) ■.^L\n,p,n) ^n(E)L^n,F), (2.1) 

Jn 

equipped with the inner product 

{ifi, = / (^(o.), V(^^))« dP(^), V^, V e 7^2. (2.2) 

We are interested in bounded operators on H that are decomposable elements 
A = (A(j)tjen, in the sense that they commute with the diagonal ones. Measur- 
able operators are defined as decomposable operators such that for all measurable 
vector's field {(p{uj),uj G ft}, the set {A{ijj)(p{u!),u! G ft} is measurable too. Mea- 
surable decomposable operators are called essentially bounded if — > is 
a element of L°°(f2,P). We set, for such A's, 



WMcm = WMoo = ess - supol|A(w)||, 



(2.3) 
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and define the following von Neumann algebra 

/C= L°°in,P, C{n)) {A: n ^ C{'H),measurable \\A\\ao < oo}. (2.4) 
Tliere exists an isometric isomorphism betwen JC and decomposable operators on 

We shall work with observables of /C that satisfy the so-called covariant property. 
Definition 2.1. A E K. is covariant iff 

UaA{uj)Ul = A{TaUj) ,yae Z,\Juj (2.5) 

We set 

^oo — {A (z IC,A is covariant}. (2-6) 

If Ua ■— Ua<E)T{~a),with the slight notation abuse where we note r for the action 
induct by r on and U = {Ua)aez, we note that 

/Coo - {Ae IC, Va e Z, [A, Ua] = 0} (2.7) 
= icn{uy, (2.8) 

so that ICao is again a von Neumann algebra. 

2.2. Noncommutative integration. The von Neumann algebra /Coo is equipped 
with the faithfuU, normal and semi-finite trace 

r{A):=E{tT{xeA{oj)xe)}, (2.9) 

where "tr" denotes the usual trace on the Hilbert space H.In the usual context of 
the Anderson model this is nothing but the trace per unit volume, by the Birkhoff 
Ergodic Theorem, whenever 7~(|v4|) < oo, one has 

r(A)= lim -Ltr(xA,A„XAj, (2.10) 

|At|->oo \Al\ 

where Al C Z and XA^ — X^aeA^ ^a- There is a natural topology associated to 
von Neumann algebras equipped with such a trace. It is defined by the following 
basis of neighborhoods: 

Nie,S) ^{Ae /Coo, 3P G /C'^™-'', \\AP\\^ < e ,T{P^) < S}, (2.11) 

where /C^°^ denotes the set of projectors of /Coo- It is a well known fact that 

A e N{e, S) ^ r(x]e,oo[(|A|)) < S. (2.12) 

As pointed out to us by V. Georgescu, this topology can also be generated by the 
following Frechet -norm on K^oo [Geo] '. 

\\A\\r= inf max{pP||oo,r(P^)}. (2.13) 

Let us denote by A^(/Coo) the set of all T-measurable operators, namely the com- 
pletion of /Coo with respect to this topology. It is a well established fact from 
noncommutative integration that 

Theorem 2.2. A^(/Coo) is a Hausdorff topological *-algebra. , in the sens that all 
the algebraic operations are continuous for the r-measure topology. 

Definition 2.3. A linear subspace £ ^ % is called T-dense if V(5 G M+, there 
exists a projection P G /Coo such that PH C £ and T{P^) < S. 
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It turns out that any element A of A1(/Coo) can be represented as an (possibly 
unbounded) operator, that we shall still denote by A, acting on H with a domain 
D{A) — {(fi G T-L, Aip G T-L] that is T-dcnsily defined. Then, adjoints, sums and 
products of elements of M.{lCoo) are defined as usual adjoints, sums and products 
of unbounded operators. 

For any < p < oo, we set 

(/Coo) := {x G K^, T{\x\P) < oo}"'"'' = {x e M{Koo), < oo}, (2.14) 

where the second equality is actually a theorem. For p > 1, the spaces Lp(/Coo) are 
Banach spaces in which LP'°(/Coo) := L*'(/Coo) n /Coo are dense by definition. For 
p = oo, in analogy with the commutative case, we set L°°(/Coo) = /Coo- 

Noncommutativc Holder inequalities hold: for any <,p, (7, r < oo so that + 
q-i = f-i, if >1 G LP(/Coo) and B G L9(/Coo), then the product AB G M{K:^) 
belongs to l7{Kao) with 

\\AB\\r<\\AUB\\,. (2.15) 

In particular, for all A G L°°(/Coo) and B G L^'(/Coo), 

\\AB\\p < P||oo||B||p and \\BA\\p < \\A\\oo\\B\\p , (2.16) 

so that L''(/Coo)-spaces are /Coo two-sided submodules of At (/Coo)- As another con- 
sequence, bilinear forms L^''°(/Coo) xL5'°(/Coo) 9 {A, B) i-> T{AB) G C continuously 
extends to bilinear maps defined on L^'(/Coo) x L^(/Coo)- 

Lemma 2.4. Let A G L^'(/Coo), P & [1, oo[ be given, and suppose T{AB) = for 

all B G L«(/Coo); P'^ + cr^ = 1. Then A = Q. 

The cascp = 2 is of particular interest since L-^(/Coo) equipped with the sesquilin- 
ear form {A, B)i^2 = T{A*B) is a Hilbert space. The corresponding norm reads 

\\A\\l = f tr(xeA:A.Xe)dP(w) = / P„Xe||^dP(a;). (2.17) 

(Where || • II2 denotes the Hilbert-Schmidt norm.) From the case p = 2, we can 
derive the centrally of the trace. Indeed, by covariance and using the centrality of 
the usual trace, it follows that T{AB) = T{BA) whenever A,Bg /Coo- By density 
we get the following lemma. 

Lemma 2.5. Let A G Lp(/Coo) and B G L«(/Coo), P~^ + = 1 be given. Then 

T{AB) = T{BA). 

Wc shall also make use of the following observation. 

Lemma 2.6. Let A G LP(/Coo) and (i?„) a sequence of elements of that con- 
verges strongly to B & /Coo- Then AB^ converges to AB in Lp(/Coo)- 

2.3. Commutators of measurable covariant operators. Let if be a decom- 
posable (unbounded) operator affiliated to /Coo with domain 2?, and A G A^(/Coo)- 
In particular H need not be T-measurable, i.e. in A^(/Coo)- If there exists a T- 
dense domain V such that AV c V, then HA is well defined, and if in addition 
the product is T- measurable then we write HA G A^(ACoo)- Similarly, if V is T- 
dense and the range of HV C 2?(A), then AH is well defined, and if in addition 
the product is T-measurable then we write AH G A^(/Coo)- 

Remark 2.7. We define the following (generalized) commutators: 
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(i) : IfAeM (/Coo) and B £ /Coo, 

[A,B]:= AB - BAe M{1C^), [B , A] := -[A, B]. (2.18) 

(ii) : If A^ LP(/Coo), B e L'^{ICoo), P,q>l such that p-^ + q-^ ^ 1 , then 

[A,B]:= AB - BAeh\]C^). (2.19) 

Definition 2.8. Let HrilCoo(i-e H affiliated to ICoo)- If A ^ A^(/Coo) *s such that 
HA and AH are in Ai{ICoo), then 

[H,A]:^ HA- AH e M{lCoo) ■ (2.20) 

We shall need the following observations. 

Lemma 2.9. 1 ) For any A e l7{JCoo),B e L«(/Coc), p,q>l, p^^ + q^^ = I, and 
Cuj e /Coo 7 we have 

T{[C,A]B} = T{C[A,B]}. (2.21) 

2) For any A, B ^ /Coo md C G L^(/Coo); we have 

T{A[B,C]] = r{[A,B]C]. (2.22) 

3) Letp,q>\ he suchthatp^^+q^^ ^ \. For any A eL^ilCoo), resp. B eL'^{ICoo), 
such that [H,A] e L^'(/Coo), resp. [H,B] £ L''{ICoo), we have 

r{[H,A]B} = -T{A[H,B]}. (2.23) 

2.4. Difierentiation. A ^-derivation 3 is a *-map defined on a dense sub-algreba 
of /Coo and such that: 

• diAB) = d{A)B + Ad{B) 

• d{A + XB) = d{A) + XdiB) 

• dlA*) = d{A)* 

• [aa, 9] = in the sense that aa o d{A) —do aa{A) Va E Z ,\/A £ /Coo- 

If di,...,dd are ^-derivations we define a non-commutative gradient by V := 
(9i, ...,9d), densily defined on /Coo- We define a non-commutative Sobolev space 

W^'P{ICoo) ■■= {A £ LP(/Coo), VA £ LP(/Coo)}- (2-24) 

and a second space for HrjK-oo, 

2?W(^) = MeLf(/Coo), HA,AH eLP{IC^)}. (2.25) 

3. The setting: Schrodinger operators and dynamics 

In this section we describe our background operators and recall from jBoGKS] the 
main properties we shall need in order to establish the Kubo formula, but within the 
framework of noncommutative integration when relevant (i.e. in Subsection 13.21) . 

3.1. Magnetic Schrodinger operators and time-dependent operators. 

Throughout this paper we shall consider Schrodinger operators of general form 

H = H{A,V) = {-iV - Af + V on L'^{R'^), (3.1) 

where the magnetic potential A and the electric potential V satisfy the Leinfelder- 
Simader conditions: 

• A(x) £ Lf„,(M'^;M'') with V - A{x) £ L^^cCI^'')- 
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• V{x) = V+{x) - V-^{x) with V±{x) e L2,c(IR'*)' ^±i^) > 0' and V-{x) 
relatively bounded with respect to A with relative bound < 1, i.e., there 
are < a < 1 and /? > such that 

\\V^ij\\ < aWAiPW + PW^PW for all e 2?(A). (3.2) 

Leinfelder and Simader have shown that H{A, V) is essentially self-adjoint on 
Cc°°(K'') [LSI Theorem 3]. It has been checked in jBoGKSj that under these hy- 
potheses H{A, V) is bounded from below: 

H(A, V)>- =: -7 + 1, so that iJ + 7 > 1. (3.3) 

(1 - a) 

We denote by xj the multiplication operator in L^(R'^) by the j"^ coordinate 
ccj, and x (xi, • • • x^). We want to implement the adiabatic switching of a 
time dependent spatially uniform electric field E^(t) • x = e''*E(t) • x between time 
t = —00 and time t = to. Here > is the adiabatic parameter and we assume 
that 

/ " e''*|E(t)|dt < 00. (3.4) 
^ — 00 

To do so we consider the time-dependent magnetic potential A{t) = A + F,f{t), 
with Fjj{t) — J^^Fjj{s)ds. In other terms, the dynamics is generated by the 
time-dependent magnetic operator 

H{t) = (-zV - A - F^(t))2 + V{x) - H{A{t), V) , (3.5) 

which is essentially self-adjoint on C^(M'*) with domain T) :— 'D{H) — V{H{t)) for 
all t e R. One has (see |BoGKSl Proposition 2.5]) 

H{t)=H -2F^{t)-D{A)+Fr,{tf onV{H), (3.6) 

where D — D(A) is the closure of (— iV — A) as an operator from L^(]R'^) to 
L2(M'^;C'') with domain C;?°(R'^). Each of its components ==Dj(A) = (-igf-- 

Aj), j — 1, . . . , d, is essentially self- adjoint on C^(K'^). 

To see that such a family of operators generates the dynamics a quantum particle 
in the presence of the time dependent spatially uniform electric field E^(t) ■ x, 
consider the gauge transformation 

[G(0^](a;) :=e'^''(*)-"i^(x) , (3.7) 

so that 

H{t) = G{t)[{-iV - Af + V]G{t)* . (3.8) 
Then if ip(t) obeys Schrodinger equation 

idtm = Hit)ij{t), (3-9) 

one has, formally, 

idtG{ty^{t) = [(-iV -Af + V + E„(t) • x] G{t)*il;{t). (3.10) 
To summarize the action of the gauge transformation we recall the 
Lemma 3.1. [BoGKSl Lemma 2.6] Let G{t) he as in dXT]). Then 

G(t)V ^ V, (3.11) 
H{t) = G{t)HG{t)\ (3.12) 
D(A + F^(i)) = D(A)-F,;(i) = G(t)D(A)G(t)*. (3.13) 
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Moreover, i[H{t),Xj] = 2D(A + F^(t)) as quadratic forms on V n V^Xj), j = 
l,2,...,d. 

The key observation is that the general theory of propagators with a time depen- 
dent generator jYj Theorem XIV. 3.1] applies to H{t). It thus yields the existence 
of a two parameters family U (t, s) of unitary operators, jointly strongly continuous 



in t and s, that solves the Schrodinger equation. 

U{t,r)U{r,s) = U{t,s) (3.14) 

U{t,t) = / (3.15) 

U{t,s)V = V, (3.16) 

idtU{t, s)ij = H{t)U{t, s)il) for ah V G 2? , (3.17) 

idsU{t, s)^ = - U{t, s)H{s)il) for aXltl^eV. (3.18) 

We refer to [BoGKSl Theorem 2.7] for other relevant properties. 



To compute the linear response, we shall make use of the following "Duhaniel 
formula". Let [/("'(t) = e~"^. For aW ij; (z V and i, s £ M we have |BoGKSl 
Lemma 2.8] 

U{t,s)%lj ^U^^\t-s)i}+i I C/(")(t-r)(2F^(r)-D(A)-F,,(r)2)C/(r,s)V'dr. (3.19) 

J s 

Moreover, 

lim U(t, s) = (t - s) strongly . (3.20) 

E|^0 

3.2. Adding the randomness. Let (ri,P) be a probability space equipped with 
an ergodic group {tq; a £ Z''} of measure preserving transformations. We study 
operator-valued maps A: 3 w H- . 

Throughout the rest of this paper we shall use the material of Section [5] with 
T-L = L^(]R'') and Z = 'L'^. The projective representation of Z'^ on H. is given by 
magnetic translations ([/(a)V')(a;) — e^°-'^^ip{x — a), S being a given d x d real 
matrix. The projection Xa is the characteristic function of the unit cube of R*^ 
centered at a £ Z"*. 

In our context natural ^-derivations arise: 

djA :^i[xj, A], j = !,■■■ ,d, Vyl = i[x,A]. (3.21) 

We shall now recall the material from |BoGKS[ Section 4.3]. Proofs of assertions 
are extensions of the arguments of |BoGKS) to the setting of LP(/Coo)-spaces. We 
refer to |Do] for details. 

We state the technical assumptions on our Hamiltonian of reference H^^ . 

Assumption 3.2. The ergodic Hamiltonian uj H> is a measurable map from 
the probability space (f2,P) to self-adjoint operators on % such that 

= H{A^,V^) = (-zV - A^f + , (3.22) 

almost surely, where A^^ (Kj j o,re vector (scalar) potential valued random variables 
which satisfy the Leinf elder- Simader conditions (see Subsection \3.1\) almost surely. 
It is furthermore assumed that is covariant as in ()2.5p . We denote by H the 
operator (-fftj)cjgn acting on %. 
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As a consequence ||/(i?cj)|| < ||/||oo and f{H) G /Coo for every bounded Borel 
function / on the real line. In particular H is affiliated to /Coo- For P-a.e. ui, 
let Ui^{t,s) be the unitary propagator associated to H^^ and described in Subsec- 
tion |33] Note that {Uuj{t, s))uj)En G /Coo (measurability in uj follows by construction 
of U^{t,s), see |BoGKS| ). For A e M{ICoo) decomposable, let 

U{t,s){A):^ r Ujt,s)A^U^{s,t)dn<-^). (3.23) 
Jn 

Then U{t, s) extends a linear operator on A^(/Coo)j leaving invariant A^(/Coo) and 
Lf(/Coo),pe]0,oo], with 

U{t,r)U{r,s) = U{t,s), (3.24) 

U{t,t) = /, (3.25) 

{U{t,s){A)Y = U{t,s){A*). (3.26) 

Moreover, lA{t, s) is to a unitary on L^(/Coo) and an isometry in LP(/Coo), P G [1, oo]. 
In addition, U(t, s) is jointly strongly continuous in t and s on LP(/Coo), P G 1, oo]- 
Pick p > 0. Let A e LP(/Coo) be such that H{ro)A and AH{ro) are in LP(/Coo) 
for some rg G [—00,00). Then both maps r H> U{t,r){A) G LP(/Coo) and i H> 
Z^(<,r)(yl) G LP(/Coo) are differentiable in LP(/Coo), with (recalling Definition [23) 



idrU{t,r){A) ^ ^U{t,r){[H{r),A]). (3.27) 
idtU{t,r){A) = [i?(i),iY(t,r)(A)]. (3.28) 

Moreover, for to < (x given, there exists C = C{to) < 00 such that for all t,r < to, 
II {H{t)+^)ll{t,r){A)\\, < C\\{H{r)+^)A\\,, (3.29) 
\\[H{t)Mt,r){A)]\\ <Ci\\{H{r) + j)A\\, + \\A{H{r) + ^)\\,) . (3.30) 



We note that in order to apply the above formula and in particular p. 271) and 
p.28p . it is actually enough to check that {H{ro) + 'y)A and A{H{ro) + 7) are in 
LJ'(/Coo). 

Whenever we want to keep track of the dependence of Uuj{t,s) on the electric 
field E — E,,(t), we shall write t/(j(E, t, s). When E = 0, note that 

i7„(E = 0,t,s) = C/i,°)(t-s) :=e-*(*-")^-. (3.31) 
For A G A^(/Coo) decomposable, we let 

U^'Hr){A):^ rC/io)(r)A^[/(°)(-r)dP(c.). (3.32) 
Jn 

We still denote by U''°^r){A) its extension to M{ICoo)- 

Proposition 3.3. Letp > 1 be given. L(^^\t) is a one-parameter group of operators 
on A^(/Coo), leaving LP(/Coo) invariant. h{^^\t) is an isometry on LP(/Coo), and 
unitary if p — 2. It is strongly continuous on LP(/Coo). We further denote by Cp 
the infinitesimal generator oflA^'^\t) in LP(/Coo).' 

{t) = e-"^" for allteR. (3.33) 

The operator Cp is usually called the Liouvillian. Let 

2?(o) = {A G L?'(/Coo), HA,AH eLP{ICoo)}. (3.34) 
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Then V^p ^ is an operator core for Cp (note that C2 is essentially self-adjoint on 

Cp{A) = [H, A] for all A e vf\ (3.35) 

Moreover, for every G fCoo there exists a sequence Bn^^ £ such that Bn^u ^ 
Bi^ as a bounded and F-a.e. -strong limit. 

We finish this hst of properties with the following lemma about the Gauge trans- 
formations in spaces of measurable operators. The map 

g{t){A) = Git)AG{ty , (3.36) 

with G{t) e'-'^-^^ ^"(^'-^ as in (|3Jl) . is an isometry on LP(/Coo), for p g]0, 00]. 

Lemma 3.4. For any p g]0,cxd], the map Q(t) is strongly continuous on h'^ICoo), 
and 

lim Q{t) = I strongly (3.37) 

>— 00 

on L^'(/Coo)- Moreover, if A £ W"^'''(/Coc)j then Q{t){A) is continuously differen- 
tiable in LP(A^oo) with 

dtg{t){A) = E„{t) ■ \7{g{t){A)). (3.38) 

4. Linear response theory and Kubo formula 

4.1. Adiabatic switching of the electric field. We now fix an initial equilibrium 
state of the system, i.e., we specify a density matrix which is in equilibrium, so 
[HulXu:] = 0. For physical applications, we would generally take Clj = f{Huj) with 
/ the Fermi-Dirac distribution at inverse temperature /3 £ (0, 00] and Fermi energy 
Ef e M, i.e., f{E) = ^^^,(Lf,) if /3 < 00 and fiE) - X(-oo,e,]{E) if /? = 00; 
explicitly 



- J — 1+o<3(Hc.-Ef) ' P^O", 

^^uj •— X(-oa,EF\\tiui) 1 P 



However we note that our analysis allows for fairly general functions / |BoGKS) . 
We set ( = (C(.j)cjGn & ^00 but shall also write instead of (. That / is the Fermi- 
Dirac distribution plays no role in the derivation of the linear response. However 
computing the Hall conductivity itself (once the linear response performed) we shall 

( E ) 

restrict our attention to the zero temperature case with the Fermi projection P,^ . 

The system is described by the ergodic time dependent Hamiltonian H^{t), as 
in (j3.5p . Assuming the system was in equilibrium at t = —od with the density 
matrix Ptj(— 00) = (^ui, the time dependent density matrix g^iit) is the solution of 
the Cauchy problem for the Liouville equation. Since we shall solve the evolution 
equation in LP(/Coo), we work with H{t) = (H;^(f))i^esi, as in Assumption 13.21 

The electric field Ej,(f) • x = e''*E(t) • x is swichted on adiabatically between 
t — —00 and t — to (typically to — 0). Depending on which conductivity on is 
interested, one may consider different forms for E(i). In particular E(t) = E leads 
to the direct conductivity, while E(t) = cos{iyt)Ei leads to the AC-conductivity at 
frequency 43. The first one is relevant for studying the Quantum Hall effect (see 
subsection 14. 4p . while the second enters the Mott's formula |KLP1 |KLM| . 



^The AC-conductivity may be better defined using the from 114.2311 as argued in |KLM| . 
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We write 

at)^G{t)CGitr ^gm), i.e., cw-/(^w)- (4.2) 

Theorem 4.1. Let 77 > and assume that e''''|E(r)|dr < 00 for all t G M. Let 
p G [l,oo[. Assume that C G VV"'^'^(A^oo) and that G V". The Cauchy problem 

( idtg{t) = [H{t),gm 

\ limt^_oo git) = C ' ^ ' 

has a unique solution in LP(A^oo), that is given by 

git) - lim Uit,s)iC) (4.4) 

s— ^ — 00 

= lim Uit,s)iCis)) (4.5) 

s— ^ — 00 

= Cit)- f dre'"'W(i,r)(E(r) • VC(r)). (4.6) 



^.(^)=^^(^,s)(g(5)), ||^.(t)||p = ||c||p, (4.7) 

for all t, s. Furthermore, git) is non-negative, and if ( is a projection, then so is 
git) for all t. 

Remark 4.2. // the initial state C is of the form (14. ip . then the hypotheses of 



Theorem \4-l\ hold for any p > 0, provided C.^ — pL^''^ that Ep lies in a region 



of localization. This is true for suitable A(^,VLj and Ep, by the methods of, for 
example, pll IWl iGKll [GK2l [GK8l IBoGKI lAENSSI iTJl ICxrHK] and for the models 
studied therein as well as in fCHj IGK3) . The bound E[||x|C(^xo||^ < oo or equiv- 
alently G L^(/Coo) is actually sufficient for our applications. For p — 1,2, we 
refer to [BoGKSl Proposition 4.2] and |BoGKSl Lemma 5.4] for the derivation of 
these hypotheses from known results. 



Proof of Theorem \41\ Let us first define 

git,s):^Uit,s)iCis)). (4.8) 
We get, as operators in M{ICoo), 

dsgit,s) = iUit,s)i[H{s),Cis)])+Uit,s)ii:^is)-VCis)) 

= ll{t,s){E„{s)-VC{s)) , (4.9) 
where we used (|3.27p and Lemma [331 As a consequence, with E,,(r) = e''''E(r), 

git,t) - git,s) ^ f dre'"'Z^(t,r)(E(r) • VC(r)). (4.10) 



Since r)(E(r) • V(C(r))||p < Cd|E(r)| ||VC||p < oo, the integral is absolutely 

convergent by hypothesis on E^(t), and the limit as s — — oo can be performed in 
LP (/Coo). It yields the equality between (H3)) and (ITO . Equality of and (H3)) 
follows from Lemma 13.41 which gives 

C= lim C(s) inLP(/Coo). (4.11) 

Since U{t,s) are isometries on LP(/Coo), it follows from (|4.4p that ||^?(i)||p = ||Cllp. 
We also get g{t) = g{t)* . Moreover, (|4.4p with the limit in lJ'{K.oo) implies that 
git) is nonnegative. 
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Furthermore, if C = then g{t) can be seen to be a projection as follows. Note 
that convergence in implies convergence in A4{JCoo), so that, 



g{t) = lim(^) U{t, s) (C) = lim^^) U{t, s) (C) U{t, s) (C) 

lim(^) U{t, s) (C) 1 ( lim(^) U{t, s) (C) 



g{t)\ (4.12) 



where we note lim'^' the limit in the topological algebra Ai{ICoo)- 

To see that g{t) is a solution of (|4.3p in LP{ICoo), we differentiate the expression 
(gSl) using (13:^8)) and Lemma lOl We get 



idtg{t) 



dre'"'W(i,r) (E(r) • VC(r)) 



dre'"'Z^(t,r) (E(r) • VC(r)) 



(4.13) 
(4.14) 



= H{t)Ac{t)- 

L — OO J _ 

= [i^(^),^'W]• (4-15) 

The integral in (|4.13p converges since by (\3.30\i , 

II [i/(t),^(t,r) (E(r) • VC(r))] ||p < 2C|| (ff + 7)(E(r) • VC)||p. (4.16) 

Then we justify going from ()4.13p to ()4.14p by inserting a resolvent {H{t) +7)^^ 
and making use of p.29p . 

It remains to show that the solution of (14. 3p is unique in L^'(/Coo)- It suffices 
to show that if i'{t) is a solution of (|4.3p with C then ^{t) = for all t. We 
define £'(^)(i) = U{s,t){v{t)) and proceed by duality. Since p > 1, with pick q s.t. 
+ (7^^ = 1. If A e 2?g°'', we have, using Lemma [2.51 

idtT[Av^'\t)]=idtT{U{t,s){A)v{t)} 

^r{[H{t),U{t.s){A)] v{t)} + r {U{t,s){A)C,{t){u{t))} 
= -r{U{t, s){A)C,{t){um + T{U{t, s){A) C,{t){um = . 

We conclude that for all t and A G Pg"-* we have 



(4.17) 



T 



[Av^'\t)^ ^t[ap^'\s)^ = T{Av{s)} . 



(4.18) 



Thus i>(^)(t) = v{s) by LemmaEll that is, v{t) ^U{t,s){v{s)). Since by hypothesis 
lim^^-oo v{s) = 0, we obtain that 1/(1) = for all t. □ 

4.2. The current and the conductivity. The velocity operator v is defined as 

v = v(A) =2D(A), (4.19) 

where D = D(A) is defined below (IX^ . Recall that v = 2(-iV - A) = i[H,x] on 
C^{Rdy We also set D(t) = D(A + F^(t)) as in (jXHl) . and v(t) = 2D(t). 

From now on g(t) will denote the unique solution to (j4.3p . given explicitly in 
gSl). If H[t)g(t) e LP (/Coo) then clearly 'Dj{t)g{t) can be defined as well by 

T>,{t)g{t) = {T»,{t){H{t) + 7)-') {Wt) + ^)g{t)) , (4.20) 

since T>.j{t){H{t) + 7)"^ e /Coo, and thus T>j{t)g{t) e LP(/Coo). 
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Definition 4.3. Starting with a system in equilibrium in state C,, the net current 
(per unit volume), J(?7, E;^,to) G , generated by switching on an electric field E 
adiabatically at rate r/ > between time —oo and time to, is defined as 

J(77, E; C, to) = r {^^{to)gito)) - r K) . (4.21) 

As it is well known, the current is null at equilibrium: 

Lemma 4.4. One has T(DjC) = for all j = 1, ■ ■ ■ ,d, and thus T (vC) = 0. 

Throughout the rest of this paper, we shall assume that the electric field has the 
form 

E(i) = f (i)E, (4.22) 

where E G gives the intensity of the electric in each direction while \£{t)\ — 0{1) 
modulates this intensity as time varies. As pointed out above, the two cases of 
particular interest are £{t) = 1 and £{t) = cos{i't). We may however, as in jKLM) . 
use the more general form 

S{t) = / cos{iyt)i{iy)di^, (4.23) 
Jr 

for suitable £{iy) (see [KLMp . 

It is useful to rewrite the current (I4.2ip . using (|4.6I) and Lemma [4.41 as 

J(77,E;C,to) = r{2D(0)(e(to)-C(io))} (4.24) 
= -r(2 r dre'"^D(0)Z^(io,r)(E(r).VC(r)) 

(-to 



= -r\2 [" dre'"^£:(r)D(0)W(to,r)(E- VC(r)) 



The conductivity tensor a{ri;(^,to) is defined as the derivative of the function 
J{t], E; C, io) : R'' ^ K'' at E = 0. Note that a{rj; (, to) is a dxd matrix {ajk (r?; (, to)}- 

Definition 4.5. For > and to G K, the conductivity tensor a{rj; to) is defined 
as 

a{r,;C,to) = M3{v,'E;C,to))\E=o , (4-25) 
if it exists. The conductivity tensor a {(^^ to) is defined by 



a(C,to) :=lima(77;C,to), (4.26) 



whenever the limit exists. 



4.3. Computing the linear response: a Kubo formula for the conductiv- 
ity. The next theorem gives a "Kubo formula" for the conductivity at positive 
adiabatic parameter. 

Theorem 4.6. Let rj > 0. Under the hypotheses of Theorem \4.1\ for p = 1, the 
current 3{r],'E; (^,to) is differ entiable with respect to E at E = and the derivative 
a(r]-X) is given by 



^jkimUo) = -r(2 r° dre'"-£(r)D,Z^(o)(to-r)(9fc(C)) 

L J —oo 

The analogue of [BESl Eq. (41)] and [SB2l Theorem 1] then holds: 



(4.27) 
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Corollary 4.7. Assume that £[t) = Ke'"*, e M, then the conductivity ajkirj; C; i^) 
at frequency v is given by 

a,k{n;C;i';0) = -r{2D, (^A +77 + zj.)-i (9fcC)}} , (4.28) 

Proof of corollary \4.7\ Recall (4.11) , in particular ( = C^C^- It follows that 
cr(r7; I/; ^; 0) in (4.27) is real (for arbitrary C, = f{H) write f — f+ — /-)• As a 
consequence , 

a(7y;j.;C;0) = -5Rr(2 T dre''V^'-D, (to - r) (9^(0)1 . (4.29) 



Integrating over r yields the result. □ 

Proof of Theorem \4 (i\ For clarity, in this proof we display the argument E in all 
functions which depend on E. From (|4.24p and Jj(?7, 0; C, ^o) = (Lemma 14. 4p . we 
have 

a,fe(77;C,io) = -|m 2r|y " dr e'"-£:(r)DEj (O)W(E, 0, r) (^^((E, r)) | . (4.30) 

First understand we can interchange integration and the limit E — > 0, and get 

^jkimCM) = -2 /" dre'"^f(r) hm r{D,(E,0)^Y(E,0,r)(afcC(E,r))} . (4.31) 

The latter can easily be seen by inserting a resolvent {H{t) +7)"^ and making 
use of dSSni), the fact that iJVC G l2{K.oo)-, the inequality : \T{A)\ < T{\A\) and 
dominated convergence. 

Next, we note that for any s we have 

lim e(E, s) = I strongly in V^{1C^) , (4.32) 

E— ^0 



which can be proven by a argument similar to the one used to prove Lemma 13.41 
Along the same lines, for B G /Coo we have 

lim g{-E,s){B^) = strongly in U, with ||e?(E, s)(B)||oo = ||B||oo • (4.33) 

-E->0 

Recalling that Dj,^(E,0) = Dj,^ - Fj(0) and that \\dkC{^,r)\\i = ||9/cClli < 00, 
using Lemma [2.61 

limr{Dj(E,0)Z^(E,0,r)(9feC(E,r))}= lim T {DjC/(E, 0, r)(afeC)C/(E, r, 0)} 

E — >-0 E — >-0 

= lim r|D,C/(E,0,r)(9fcC)C^^°H'')| , (4-34) 

where we have inserted (and removed) the resolvents (i?(E, r)+^)^^ and (i/+7)^^. 

To proceed it is convenient to introduce a cutoff so that we can deal with Dj 
as if it were in /Coo- Thus we pick G C^(M), real valued, |/„| < 1, /„ — 1 on 
[— n,n], so that fn{H) converges strongly to 1. Using Lemma [2.61 we have 

r{D,C/(E,0,r)(9feC)C/^°H^)} = „l™^r{/„(i?)D,C/(E,0,r)(9,C)C/^°HO} 

= lim T{u{E,0,r){{dkCKH + j))U^°HrKH + jr'fn{H)-DA 

= T{u{E,0,r)mOiH + 7)){u^'Hr)(H + ^)-'B,)} , (4.35) 
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where we used the centrahty of the trace, the fact that {H + 7)^^ commutes with 
?7^°'' and then that {H + j)~^'Dj G /Coo in order to remove to hmit n ^ 00. Finahy, 
combining (j4.34p and (|4.35l) . we get 

Um r {Dj (E, 0)U(E, 0, r) (<9fcC(E, r))} (4.36) 

E — >0 

= r{D,W(°)(-r)(9feC)}. (4.37) 

The Kubo formula (|i:?fl) now follows from (|i3T|) and (|T77)) . □ 

4.4. The Kubo-Stfeda formula for the Hall conductivity. Following |BES[ 
lAGj . we now recover the well-known Kubo-Stfeda formula for the Hall conductivity 
at zero temperature (see however Remark 14.111 for AC-conductivity). To that aim 
we consider the case £{t) = 1 and Iq — 0. Recall Definition 14.51 We write 



= ^.-.^(^^''"^0) , and af/'(r;) = a,- ,(^; 0) . (4.38) 

Theorem 4.8. Take £(t) — 1 and = 0. If C = P^^^-* is a Fermi projection 
satisfying the hypotheses of Theorem \4-l \ with p ~2, we have 

ij-^p{EF) djP^'^'-\dkP^^'''> }, (4.39) 



for all j, fc = 1, 2, . . . , d. As a consequence, the conductivity tensor is antisymmetric; 

J, 3 



in particular a) = for j = 1,2, . . . ,d. 



Clearly the direct conductivity vanishes, cj^'^-' — 0. Note that, if the system is 
time-reversible the off diagonal elements are zero in the region of localization, as 
expected. 



Corollary 4.9. Under the assumptions of Theorem \4-S\ if A. ~ Q (no magnetic 
field), we have cj^'^^ = for all j,k = 1,2, . . . , d. 

We have the crucial following lemma for computing the Kubo-Stfeda formula, 
which already appears in jBESj (and then in |AGj ) . 

Lemma 4.10. Let P G IC°° be a projection such that dkP G \J'{1C°°), then as 
operators in M.{K.°") (and thus in U'{JC°°)), 

dkP=[P,[P,dkP]]. (4.40) 

Proof. Note that d^P = dkP^ = PdkP + {dkP)P so that multiplying left and right 
both sides by P imphes that P{dkP)P = 0. We then have, in LP(/C°°), 

dkP - PdkP+idkP)P = PdkP+{dkP)P-2P{dkP)P 
= P{dkP){l - P) + {1 - P){dkP)P 
= [P,[P,dkP]]. 

□ 

Remark that Lemma (|4.10p heavily relies on the fact P is a projection. We shall 
apply it to the situation of zero temperature, i.e. when the initial density matrix is 
the orthogonal projection p(-^^). The argument would not go through at positive 
temperature. 
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Proof of Theorem \4-8\ We again regularize the velocity Dj.ui with a smooth func- 
tion /„ G C^(R), \fn\ < 1, fn — 1 on [— n,n], but this time we also require that 
/„ = outside [-n - l,n+ 1], so that /„X[-n-i,«+i] = fn- Thus 'Djfn{H) G 
LP^°ilCao), < p < oo. Moreover 

/„(i?)(2D,)/„(i/) = fn{H)Pn[2T>,)PMH) = -f^{H)dj{PnH)f^{H) (4.41) 

where P„ = P^ — X[-n-i.n+i]{H) so that i/P„ is a bounded operator. We have, 
using the centrality of the trace T , that 

?(f)(r) := -r{2D,-^Z^W(-r)(9,p(^-))} (4.42) 
= - lim r|w(°)(0(/4^)2D,,^/„(i/))afcP(^-)|. (4.43) 
Using Lemma \2M and applying Lemma [4.101 applied to P — pi^^) ^ it follows that 
T{u^°Hr)ifniH)2D,fr,iH))dkP^^-^} (4.44) 

= r{z^(°'(^)(/nW2D,/„(i/)) [p^^-), [p(^-),afcP(^-)]]} 

= r{z^(°'(r) ([p(^-), [p(^-),/„(i/)2D,/„(i/)]]) afep(^-)} , 

= -r{w(")(r) ([p(^-),/„(77) [p(^-),5,(i/P„)] /„(i/)]) 9,P(^-)}, 

where we used that p(-^^) commutes with J/'"-* and /„(i?), and (|4.4ip . Now, as 
elements in M{IC°°), 



P^^^\djHPn = HPn.djP^^^"^ 



(4.45) 



Since [H, djP^'^"^]] is weh defined by hypothesis, fn{H) [HP^, djP^^"'^] fn(H) con- 
verges in PP to the latter as n goes to infinity. Combining (|4.43l) . (I4.44p . and (|4.45l) . 
we get after taking n ^ oo, 



p{EF) 



Next, 



H, ap(^^) 



p, p(^^),a,p(^^) 



so that, recalling Proposition 13. 3[ 



where (A,P)l2 = T{A*B). Combining gjH]), (|442]) . and (1448)) . we get 

^T^f ^(r;) = - (z(/:2 + ^r;)-i/:2 ([p(^-),a,p(^-)]) ,afeP(^-))^^ . 

It follows from the spectral theorem (applied to £2) that 

lim (£2 + i^y^ C2 = P(Kor£2)^ strongly in l^ilCao) , 

?7— >0 ^ ' 



(4.46) 
(4.47) 

(4.48) 
(4.49) 
(4.50) 



where P(Kcr£2)^ orthogonal projection onto (Ker£2)^- Moreover, as in 

[BoGKS] one can prove that 



piEF)^Q^p(E^)'] ^ (Kcr£2 



(4.51) 
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Combining (j449l) . (|450l) . (|43T|) . and Lemma [231 we get 

which is just ((4391) . □ 



Remark 4.11. // one is interested in the AC- conductivity, then the proof above is 
valid up to (14.491) . In particular, with £{t) — 5Re"^*, one obtains 

a['^^\r^)^-^U{C2 + v + iri)-^ C2(\p'~''^\djP^''-^\] ,dkP'^''-^) , . (4.52) 



L2 



The limit r] ^ can still be performed as in [KLMl Corollary 3.4] . It is the main 
achievement of jKLMj to be able to investigate the behaviour of this limit as v ^ 
in connection with Mott 's formula. 
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